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We report the specific heat measurements on the newly discovered Fe-based layered superconductor 
LaO0.9F0.1-.5Fe As with the onset transition temperature T c ss 28 K. A nonlinear magnetic field dependence 
of the electronic specific heat coefficient y(H) has been found in the low temperature limit, which is consistent 
with the prediction for a nodal superconductor. The maximum gap value A ~ 3.4±0.5 meV was derived by 
analyzing y(H ) based on the d-wave model. We also detected the electronic specific heat difference between 9 
T and T in wide temperature region, a specific heat anomaly can be clearly observed near T c . The Debye tem- 
perature of our sample was determined to be about 315.7 K. Our results suggest an unconventional mechanism 
for this new superconductor. 

PACS numbers: 74.20.Rp, 74.25.Bt, 65.40.Ba, 74.70.Dd 



In the superconducting state, the conduction electrons pair 
and condense into a low energy state leading to the forma- 
tion of a gap at the Fermi energy £> . It is this gap that pro- 
tects the superconducting condensate. Based on the general 
space group of the material, and together with the detailed 
pairing interaction between the two electrons, the supercon- 
ducting gap should have a specific symmetry for an individual 
superconductor. This gap symmetry is very important for un- 
derstanding the underlying mechanism of a superconductor. 
For example, a d x 2_ v 2 symmetry has been proved by tremen- 
dous experiments for the cuprate superconductors.] 1]. One 
effective way to detect the gap symmetry is to generate quasi- 
particles (QP) from the condensate and then trace out the de- 
tailed way for the accumulation of the QP density of states 
(DOS). Specific heat is one of the powerful tools to mea- 
sure the DOS at the Fermi level. Textbook knowledge tells 
us that, the low temperature electronic SH for an s-wave gap 
should have an exponential temperature dependence, namely 
C e oc exp(-To/T), here Tq is a characteristic temperature re- 
lated to the magnitude of the energy gap. For a supercon- 
ductor with nodal gap symmetry, the DOS is a power law of 
energy leading to a power law dependence of temperature for 
the electronic SH: C e oc T 2 for the gap with line nodes and 
C e oc r 3 for point nodes. J2] For example, in cuprate super- 
conductors, there are line nodes in the gap function, this re- 
sults in[3] an electronic SH C e = aT 2 , where a oc y n /T c and 
y„ is the specific heat coefficient reflecting the DOS at the 
Fermi level of the normal state. In the mixed state, the mag- 
netic vortices will induce depairing both within and outside 
the vortex cores leading to the localized and delocalized QP 
DOS, respectively. Volovik Q4£| pointed out that for d-wave 
superconductors in the mixed state, supercurrents around a 
vortex core lead to a Doppler shift to the quasi-particle ex- 
citation spectrum, which affects strongly the low energy ex- 
citation around the nodes. It was shown that the contribution 
from the delocalized part will prevail over the core part and 
the SH is predicted to behave as|H, 0] C vot = ky n T ^H/H c2 
with k in the order of unity, H C 2 the upper critical field. This 
prediction has been verified by many measurements which 



were taken as the evidence for d-wave symmetry in cuprate 
superconductors!! 0, 1 i M M Q ■ 

In past years, superconductivity at several Kelvins was ob- 
served in some quaternary oxypnictides with a general for- 
mula as LnOMPn (where Ln= La and Pr; M = Mn, Fe, Co 
and Ni; Pn = P and As)|d O. It was found that the 
T c can be increased by partially substituting the element O 
with F[14]. Recently, it was reported that the superconduct- 
ing transition temperature can be increased to 26 K in the ma- 
terial LalCh-jFJFeAs (x = 0.05 ~ O.12)|0. This is sur- 
prising since the iron elements normally give rise to magnetic 
moments, and in many cases they form a long range ferromag- 
netic order, and are thus detrimental to the superconductivity 
with singlet pairing. Therefore it is highly desired to know 
the gap symmetry of this new superconductor. In this Letter, 
we report the measurements on low temperature specific heat 
under different magnetic fields. A nonlinear field dependence 
of SH coefficient y is discovered. Our data together with a de- 
tailed analysis indicate that the new superconductor LaOFeAs 
may have a nodal gap. 

The polycrystalline samples were synthesized by using a 
two-step solid state reaction method. First the starting materi- 
als Fe powder (purity 99.95%) and As grains (purity 99.99%) 
were mixed in 1:1 ratio, ground and pressed into a pellet 
shape. Then it was sealed in an evacuated quartz tube and 
followed by burning at 700°C for 10 hours. Then the resul- 
tant pellet was smashed and ground together with the LaF3 
powder (purity 99.95%), La203 powder (purity 99.9%) and 
grains of La (purity 99.99%) in stoichiometry as the for- 
mula LaOo.gFo i-^FeAs. Again it was pressed into a pellet 
and sealed in an evacuated quartz tube and burned at about 
940°C for 2 hours, followed by a burning at 1150°C for 48 
hours. Then it was cooled down slowly to room tempera- 
ture. Since a little amount of F may escape during the sec- 
ond step fabrication, in the formula for our sample, we use 
(0.1 - 5) as the possible concentration of F. X-ray diffraction 
patterns in this sample show that the dominant component is 
from LaOo.yFo.i-tfFeAs. 

The ac susceptibility were measured based on an Oxford 
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FIG. 1: (color online) Raw data of specific heat plotted as C/T vs 
T 2 are shown in the main frame. The inset (a) and (b) show the 
temperature dependence of the resistivity and the ac susceptibility at 
zero dc magnetic field, respectively. 



cryogenic system (Maglab-Exa-12). The resistivity and the 
specific heat was measured on the Quantum Design instru- 
ment physical property measurement system (PPMS) with 
temperature down to 1.8 K. We employed the thermal relax- 
ation technique to perform the specific heat measurements. 
For the SH measurements, we used a latest improved SH mea- 
suring puck from Quantum Design, which has negligible field 
dependence of the sensor of the thermometer on the chip as 
well as the thermal conductance of the thermal linking wires. 

In the main frame of Fig. 1 we show the raw data of spe- 
cific heat plotted as C/T vs T 2 under two different magnetic 
fields T and 9 T. Very similar to the case in cuprate super- 
conductors, no visible superconducting jump can be seen from 
the raw data at zero field, indicating a rather small superfluid 
density or condensation energy in the present system. By ex- 
trapolating the data in the superconducting state at zero field 
down to K, it is surprisingly to see a rather small resid- 
ual term y$ « 0.69 mJ/molK 2 , indicating a very small non- 
superconducting volume or little impurity scattering centers 
in our present sample. [Here one mole means one unit cell 
or two molecules, i.e.,(LaOo.9Fo.i-5FeAs)2 ] Inset (a) of Fig. 
1 shows the temperature dependence of resistivity at T in a 
wide temperature range up to 300 K. The onset transition tem- 
perature is about 28 K. One can also see that the residual resis- 
tivity ratio (RRR) is about 18.5, which is a rather large value 
for a polycrystalline sample. The residual resistivity at 30 K 
is only 0.13 mO cm, which is several times smaller than the 
reports by other groups lflrjl . All these suggest that our sam- 
ples are much cleaner with fewer scattering centers or smaller 
non-superconducting volume. The ac susceptibility data mea- 
sured at zero dc magnetic field were shown in the inset (b) of 
Fig. 1. 

The raw data of the specific heat in various magnetic fields 
at T < 7.7 K are shown in Fig. 2(a). One can see that the 
Schottky anomaly is visible when fioH > 2 T, while it is neg- 



ligible at lower fields because the SH peak due to the Schot- 
tky anomaly still locates in the lower temperature region (e.g. 
lower than 1.8 K) when the field is low. So we just need to 
consider the correction from the Schottky anomaly at high 
fields when fitting the data to separate the electronic contri- 
bution from other contributions. As we all know, the total 
specific heat of a superconductor can be expressed as 



C(T,H) = [ ro + y(H)]T + {3T 3 + C Sch (T,H), 



(1) 



where the four terms represent the residual electronic spe- 
cific heat, the electronic contribution, the phonon contribution 
and the magnetic impurity contribution (the so-called Schot- 
tky anomaly), respectively. The two-level Schottky anomaly 
is given by nx 2 e x /(l + e x ) 2 (x = gfiBH/k B T) at nonzero fields, 
where jj B is the Bohr magneton and n is the concentration of 
paramagnetic centers. As revealed by the thick solid lines in 
Fig. 2(a), the data under the fields jUqH > 2 T can be fitted 
very well by equation (1). As for the data of lower fields, 
we just simply extrapolated them down to K linearly and 
consequently the electronic contribution and the phonon con- 
tribution were separated. 

From the treatment of the data as described above, we ob- 
tained the fitting parameters for different fields, as shown in 
Table-I. Here we take the average value f3 ~ 0.49 mJ/(mol K 4 ). 
And the residual term yo was determined to be 0.69 mJ/(mol 
K 2 ), which is an extremely small value compared with that in 
other systems. The small contribution from Schottky anomaly 
and residual term make our data analysis more easy and reli- 
able. Using the obtained value of ft and the relation &o = 
(l2n 4 k B N A Z/5p) 1/3 , where N A = 6.02 xlO 23 moF 1 is the Avo- 
gadro constant, Z = 8 is the number of atoms in one unit cell, 
we get the Debye temperature &o = 315.7 K for the present 
sample. Suppose it is in the weak electron-phonon coupling 
regime, we then use the McMillan equation to evaluate the 
electron-phonon coupling strength X e -ph via lfl7ll 
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where pt is the Coulomb pseudopotential taking about 0.10. 
Using 0£> = 315.7 K and T c - 26 K, we obtain a very large 
value for the electron-phonon coupling: A e - P h = 1.2. This 
may suggest that a simple model based on electron-phonon 
coupling and the McMillan equation cannot interpret the high 
transition temperature in the present system. 

In the following we investigate the magnetic field induced 
change of the SH coefficient y(H). After removing the 
Schotkky anomaly for each field, we are left with only the the 
first three terms in eq.(l). The results are shown in Fig. 2(b). 
It is clear that the low temperature part is quite straight for all 
fields, this allows to determine the zero temperature value of 
y{H) at different fields. As shown in Fig. 3, y(H) increases 
nonlinearly as the magnetic field increases from T to 9 T. 
In fact, the nonlinear behavior can be roughly described by a 
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FIG. 2: (color online) Temperature and magnetic field dependence of 
the specific heat in C/T vs T 2 plot in the low temperature region, (a) 
Raw data before removing the Schottky anomaly. We can see that the 
Schottky anomaly becomes visible when the field is higher than 2 T. 
The thick solid lines represent the theoretical fit (see text) containing 
all terms in eq.(l). (b) Replot of the data after the Schottky anomaly 
was subtracted. The dotted lines in (a) and solid lines in (b) are 
guiding to the eyes. 



TABLE I: Fitting parameters of equation (1). 
fi H(T) y(#)(mJ/molK 2 ) /3(mJ/molK 4 ) n(mJ/molK) g 
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simple equation y(H) = A V7/ as shown by the red solid line, 
which is actually the theoretical prediction for superconduc- 
tors with line nodes in the gap function[4]. This suggests that 
in the 10% F-doped LaOFeAs sample the gap clearly has a 
nodal structure. 

Although our data shows a relation being close to the d- 
wave prediction y(H) oc y/H, it is by no means to say that the 
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FIG. 3: (color online) Field dependence of the field-induced term 
y(H) at T = K (symbols). The solid red line is the fit to y(H) = 

aVh. 



gap in the present sample is definitely of d-wave type, since 
other type of pairing symmetry with nodes can also give rise 
to a nonlinear y(H). While we can nevertheless use the d- 
wave model to derive some important parameters. For exam- 
ple, we can obtain the maximum gap value by investigating 
the field-induced term y(H) quantitatively, as we have done 
successfully in LaSrCuO single crystals iflll [Till . The term 
y(H), which mainly arises from the Doppler shift of the quasi- 
particle excitation spectrum near the nodes induced by the su- 
percurrent flowing around vortices, has a direct relation to the 
slope of the gap at the node, v&=2Ao/hkf with Ao the d-wave 
maximum gap in the gap function A = Ao cos(20), kf the 
Fermi vector( taking k,F ~ n/a ~ 0.78 A -1 , here a = 4.03A is 
the in-plane lattice constant). We have known that the relation 
between Va and the prefactor A is given by 



4k B jlT riV mol a L 



V'A 



37f 



(3) 



where Oo is the flux quantum, n is the number of conduction 
planes per unit cell, d is the c-axis lattice constant, V mo i is 
the volume per mole, and cil = 0.465 for a triangular vortex 
latticed 12011 . The value of the prefactor A, which is about 
3.2 mJ/(mol K 2 T° ,s ), has been obtained from fitting the data in 
Fig. 3. So we can extract gap value Ao ~ 3.4 + 0.5 meV using 
the known parameters for our sample. This is a reasonable 
value with the ratio 2Ao/fcBr c ss 4.0, if we take T c = 20 K. 
This ratio is quite close to the prediction {2Aq I ksT c = 4.28) 
for the weak coupling d-wave superconductors. 

The nonlinear y{H) found here cannot be simply attributed 
to the multigap effect as observed in MgZ?2[21] and 2H - 
NbSe2MM in which the y{H) is nonlinear. The reason is 
that, even if there are multigaps in these two systems, the zero 
field SH data y(T) shows a clear flattening in low tempera- 
ture region corresponding to the weak excitation of QPs for 
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FIG. 4: (color online) Temperature dependence of y(0T) - y(9T) up 
to 30 K is shown in the main frame. One can see a clear specific 
heat anomaly near T c . The inset shows the temperature dependence 
of resistivity at 9 T from 2 K to 40 K. 



an s-wave superconductor. This is completely absent in the 
present system, as shown in Fig. 1 . One may further ascribe 
the nonlinearity found for LaOFeAs to the granular feature of 
the sample: the grains are randomly aligned within the bulk 
sample, when a field is applied, the creating rates of DOS by 
a magnetic field are different among the grains with different 
orientations. In this case, a non-linear y(H) may be observed, 
especially in a system with high anisotropy. This possibility 
cannot ruled out, but it is difficult to understand why the rela- 
tion y(H) oc V7/ is roughly satisfied here. 

In Fig. 4 we present the temperature dependence of y(QT) - 
y{9T), one can see that a shallow SH anomaly starts at about 
25 K and ramps up slowly with decreasing temperature, and 
shows a peak at about 15 K, which is just the middle resis- 
tive transition point at 9T. This broadened anomaly may be 
induced by both the broad transition at zero field (the super- 
conducting phase is still not perfectly uniform) and the very 
low superfluid density leading to a strong phase fluctuation. 
Future experiments on improved samples will fix these prob- 
lems. Although a magnetic field of 9 T is still difficult to sup- 
press the superconductivity completely, it is clear that the dif- 
ference between 9T and T does not show a flattening down 
to about 2.8 K. This is not expected by an s-wave BCS super- 
conductor with T c beyond 20 K. This fact may also corrob- 
orate our conclusion derived from the nonlinear behavior of 
y(H) at zero K that there are nodes on the gap function. 

In summary, the low temperature specific heat measure- 
ments reveal that the new superconductor LaOFeAs has a 
rather low superfluid density and condensation energy. The 
field induced extra DOS y(H) follows a nonlinear behavior 
which is roughly proportional to V//. This may suggest that 



LaOFeAs superconductor has a nodal gap and is probably an 
unconventional superconductor. 
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